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ABSTRACT 

The Fermi satellite has detected GeV emission from a number of gamma-ray 
bursts and active galactic nuclei at high redshift, z ^ 1.5. We examine the con- 
straints that the detections of gamma rays from several of these sources place on 
the contribution of population-Ill stars to the extragalactic background light. 
Emission from these primordial stars, particularly redshifted Lyman-a emission, 
can interact with gamma rays to produce electron-positron pairs and create an 
optical depth to the propagation of gamma-ray emission, and the detection of 
emission at >10 GeV can therefore constrain the production of this background. 
We consider two initial mass functions for the early stars, and use derived SEDs 
for each to put upper limits on the star-formation rate density of massive early 
stars from redshifts 6 to 10. Our limits are complementary to those set on a high 
near-IR background flux by ground-based TeV-scale observations, and show that 
current data can limit star-formation in the late stages of reionization to less 
than 0.5 M0 yr^^ Mpc^'^. Our results also show that the total background 
flux from population-Ill stars must be considerably less than that from resolved 
galaxies at wavelengths below 1.5 fim. 
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I INTRODUCTION 

The reionization o f the universe, which t ook place around 
a redshift of 10.5 l|Komatsu et al.l|201ll ). is generally be- 
lieved to be driven primarily by ionizing photons from 
early 'Population III' (pop-Ill) stars. As these stars form 
from primordial unenriched, metal-free hydrogen and he- 
lium, they undergo a formation process that is sub- 
stantially different from that of later population I and 

II stars. Simulations of the pr oduction of pop-Ill stars 



11 stars, oimulation s ot tne pr oduction ot pop-Ill stars 
dAbel. Bryan fc Norman 2000l: Bromm. Coppi fc LarsonI 



20021: iTan fc McKee|[200i: lYoshida et "al1l2006l : iNormanI 



20081 ) generally find an initial mass function (IMF) that 



is heavily biased towards high masses, 10 - 1000 Mq. A 
considerable fraction of the radiant energy from these 
stars is released at ionizing wavelengths (< 912 A), 
which allows reionization of the universe to be com- 

J leted on the times cale required by Lyman-a forest data 
Becker et al.ll200ll ). Despite their importance in cosmol- 
ogy and impact on IGM evolution, pop-Ill stars con- 
tinue to evade direct detection. Detecting the redshifted 
UV emission from pop-Ill stars is a primary goal of 
the upcoming James Webb Space Telescope, though even 
with the state of the art sensitivity of this instrument 
detecting individual metal-free star s will be challenging 
iRvdberg. Zackrisson fc ScottlbOlJ ). Searching for indi- 
rect evidence of these stars and their integrated cosmo- 
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logical impact is therefore the primary way of under- 
standing the properties of the reionization-era universe. 

Photon production from the reionization era is en- 
coded in the evolving spectral energy distribution (SED) 
of the accumulated photon background, which we observe 
locally as the extragalactic background light (EBL). Red- 
shifted UV radiation from pop-Ill stars can be expected 
to appear as a distinct component of the near-IR por- 
tion of the EBL, and the spectral details of this observ- 
able light could in principle inform an observer about 
the redshift of reionization and the nature of the sources 



responsible (jK ashlins kv et al. 2004; Cooray & Yoshida 
12004 iMadaufc Silkll26o5l : iFernandez fc Komatsull2006l' l. 
However, observations of the absolute intensity of the 
EBL in the near-IR are severely hindered by the presence 
of bright galactic foregrounds, which are produced by di- 
verse sources including stars and the interstellar medium 
(ISM) of the Milky Wa y and 'zodiacal light' f rom dust 
within our solar system (jHauser fc Dwelj 200 ih . 

It has been proposed that a high level of IR back- 
ground could be due to radiation from the first generation 
of stars. Multiple experiments, most notably the DIRBE 
experiment on the COBE satellite, have attempted to 
measure the sky brightness at near-IR wavelengths, and 
foreground subtraction analyses have been presented 

I 1 1 ' 1 ; 1 

by a number of authors dW right & Reese 2000; Wright' 



2001 



2001 



Goriian. Wright fc Charv 2000: C ambrcsv ct al. , 
Levenson. Wright fc Johnsc3l2007h . with extra- 



galactic components generally exceeding the contribu- 
tion o f resolvable galaxies. As discussed in lGilmore et al.l 
l|201l[ l (GSPDll), the discrepancies between these var- 
ious measurements are largely due to the models for 
the impact of zodiacal light assumed. The IRTS experi- 
ment on the Space Flyer Satellite found an even higher 
meas urement than DIRBE in the near-IR at 1.4 to 4 
Plvcl (|Matsumoto et alll2005l ). and linked this radiation 
excess to pop-Ill formation. 

Other findings have cast doubt on this interpre- 
tation of the high near-IR background flux arising 
from reionization-era sources. TeV-scale observations by 
ground-based atmospheric Cherenkov telescopes can con- 



far-IR jDwek fc KrennrichI 20051: 


Aharonian et al.ll2006l: 


iMazin fc Rauell2007l: lAlbert et al. 


12008!) through the od- 



tical depth introduced by pair-production interactions 
with background photons to gamma rays from extra- 
galactic sources such as gamma-ray emitting blazars. 
Constraints of this type, which assume a limit to the 
intrinsic hardness of gamma-ray spectra, have generally 
found agreement with a low value for the EBL in the 
near-IR. While these limits can be violated in some cir- 
cumstances, a EBL flux as high as that measured by the 
IRTS experiment and the higher DIRBE claims would 
almost certainly not be po ssible without the em ergence 
of new physics. Moreover, iMadau fc Sil3 (|2005l ) (here- 



after MS05) investigated the theoretical implication of 
the claim of a large IR background excess, and found 
that creating a near-IR background excess would be in- 
consistent with the observed chemical enrichment history 
of galaxies. 

In this work, we revisit the possibility of constrain- 
ing the contribution of reionization-era stars to the 
EBL using gamma-ray observations, and we then ad- 
dress the corresponding limits on high-redshift star for- 
mation. Rather than focusing on the relatively low- 
redshift blazars that have been viewed with ground-based 
gamma-ray telescopes and that have been the target of 
aforementioned EBL studies, here we explore the im- 
plications of the most distant sources observed by the 
Large-Area Telescope (LAT) on the Fermi Satellite at 
lower GeV energies. The gamma-ray optical depth for 
high-redshift blazars and gamma-ray bursts (GRBs) due 
to the EBL produced by galaxies in a variety of mod- 
els was the topic o f a recent Fermi collaboration paper 
IIAbdo et al. l l2010al V However, none of the EBL models 
discussed in this paper included explicitly a contribution 
from the high-redshift UV sources responsible for reion- 
ization, like p op-Ill stars. Here, w e will employ a simi- 
lar method to lAbdo et al.l (|2010al '). by using the highest 
energy photons seen from several high-redshift sources, 
together with the spectrum observed at lower energies, 
to put an upper limit on the background flux, and we 
treat the flux from pop-Ill stars as a separate component 
from the EBL produced by lat er galaxies. Our work is 
also re lated to the calculation of lRaue. Kneiske fc MazinI 
(|2009') (RKM09), which examined the limits possible on 
low- and zero-metallicity stars from gamma-ray data and 
placed limits on reionization-era star formation from red- 
shifts 7 to 15. This analysis was limited to considering 
the EBL bounds placed by comparatively low-redshift 
TeV blazars, rather than the high redshift GeV sources 
we consider in this work. 

The prospects of using high redshift observations 
with Fermi, in particular of distant gamma-ray bursts, 
as probes of light originating in t he reionization era 
was noted in iKashlinskvl (120051 ) and feashlinskv fc Bandl 
l|2007h . These papers recognized that the gamma-ray op- 
tical depth introduced by photons from early stars would 
create a (1 -I- z)~^ dependence in the gamma-ray mean 
free path, due to the constant comoving number den- 
sity of the resulting background photons, and that the 
highest-redshift sources would therefore be the most sen- 
sitive to the spectral signature of this radiation. High red- 
shift gamma-ray spectra, when used as a probe of back- 
ground photons, can therefore effectively isolate an EBL 
component originating in the early universe. Though the 
LAT has been less successful at detecting high-energy 
GRBs than envisioned by these authors, high-redshift 
detections of both GRBs and blazars have been made at 
energies above 10 GeV, and these are sufficient to con- 
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strain a large excess at the shorter bands of the near-IR 
wavelengths. 

In the next section, we first review the calcula- 
tions of the buildup of background photons and the phe- 
nomenon of gamma-ray attenuation through electron- 
positron pair-production, and then describe a model for 
the SED of the background component of pop-Ill stars. 
In Section [3l we present our results for the limits that 
can be obtained from observed gamma-ray sources on 
the portion of the local EBL that originates from high- 
redshift, followed by upper bounds on the total rate of 
pop-Ill star formation. Section [4] presents conclusions 
and the prospects of strengthening our limits with future 
high-redshift gamma-ray source detections. We assume 
a ACDM cosmology with par ameters consistent wit h a 
maximum likelihood results of iKomatsu et all (|201ll ) m 
this work: h = 0.702, Ot = 0.0455. 



2 METHODS 

To examine the effect that the photon populations cre- 
ated during reionization have on gamma-ray propaga- 
tion, we will develop template SEDs that describe the 
emissivity of the universe during these early times. These 
templates can be freely rescaled to simulate varying star- 
formation rate densities. They are described in full detail 
in ^TSl 

In the following sections, we will treat independently 
the hypothetical EBL component from reionization-era 
pop-Ill stars and the component produced primarily by 
later population I/II stars. We will distinguish between 
these by referring to the first as the reionization-era EBL 
('r-EBL') and the second as the post-reionization EBL 
('p-EBL'). We emphasize that this work is not intended 
to comment on the relationship between these compo- 
nents, or the details of how the transition from pop-Ill 
to population-II star formation occurs, but only to ex- 
amine some possible scenarios for the r-EBL, of which 
we presently have only indirect evidence. 



mological line element, defined as 

dl c 1 

dz 



(2) 



(1 + ^)Ho ya^TTIpnH 

for a fiat ACDM universe (|Peebledll99l ). An EBL com- 
ponent originating entirely above a redshift Zr will evolve 
only through passive redshifting &t z < Zr, with an invari- 
ant comoving photon number density. If this component 
is observed locally with a spectral energy distribution 
Jl{vo), then the proper fiux of these photons at redshift 

2 < 2r is 

Jl{v,z) = {z + lfjl{vo); vo^u/{z + l) (3) 

If Jl refers in this case to the EBL produced by 
early low metallicity stars, and Zr is the end of the era of 
preeminence of these sources, then it is easy to see that 
Jl{z) will become an increasingly large fraction of the 
total background as z approaches Zr- This is due to the 
fact that the majority of the background light emerging 
from res olved galaxies after rei o nization in recent mod- 
els (e.g. ISomerville et al.ll201ll: IPqrnmguez et aLllioTl] : 
iFranceschini. Rodighiero fc Vaccaril l200a i comes from 
redshifts considerab ly low er than Zr . In the fiducial model 
oflSo rncr villc et al] ()201ll l. for instance, ~ 75 per cent of 
p-EBL photons between 0.1 and 10 microns are emitted 
at < 2 < 2, and 96 per cent at 2 < 4. 

2.2 Gamma-ray attenuation 

Gamma-gamma scattering into electron-positron pairs 
can occur when there is sufficient energy in the center- 
of-mass frame of the two-photon system. Including the 
effect of interaction angle as measured in the cosmologi- 
cal frame, this condition is 



y/2ET_E2{l - cose) ^ 2meC^, (4) 

where Ei and E2 are the photon energies and 6 is the 
angle of incidence, which for our purposes is a random 
distribution on the unit sphere. We can rewrite Equation 
[4]to define the minimum threshold energy Eth for a back- 
ground photon to interact with a gamma ray of energy 

E-y, 



2.1 Evolution of background photons 

In general, the EBL in place at a given redshift zq can 
be calculated by an integral over source terms at all red- 
shifts 2 > 20, with additional terms describing the effect 
of redshift on photon wavelength and flux redshift (e.g., 
ISomerville et al.ll201ll ). 



Ji^(l^O, zo) 



— ^ —^(v,z)dz, (1) 



AtvJ^^ dz (l + z) 



where e{i>, 2) is the galaxy emissivity at redshift 2 and 
frequency 1/ = voil + z) / {1 + zq) , and dl/dz is the cos- 



Eth = 



(5) 



E-,{1 - cos 61) ■ 

The cross-section for this process i s jBreit fc Wheeleii 
1 19341 : [Gould & Schroder 1967; Madau fc Phinnevlll996r 



a{EuE2,e) = ^il-l3'') 



(6) 



2/3(/3' - 2) + (3 - In ' ^ + ^ 



where 



2m2c'' 



£;i£;2(i-cose)' 



(7) 
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and or is the Thomson scattering cross section. We can 
also revisit Eq. [5] to define the minimum energy or max- 
imum wavelength for the background photons that can 
interact with a gamma ray of energy as 



Ebg -- 

or equivalently, 



A 



- 475 



10 GeV 



A. 



(8) 



(9) 



These equations all refer to rest-frame energies. 

To calculate the optical depth for a gamma 
ray observed at energy Ery, we perform the integral 
along the line of sight to the t arget at redshift z 
l|Stecker. de Jager fc Salamonlliggd ). 

1 1""° dl 

t{E^,zo) = -I dz—j du{l-u) (10) 



X ; dEbgn{Etg,z) a{E-,(l + z),Etg, 



Where we have 

Emin ~ Eth (1-1-2:) 



2m? 



E^{1 + z){l-cose) 



to account for the redshifting of the gamma-ray energy. 
Here n{Ebg, z) is the proper density of target background 
photons as a function of energy Etg and redshift z, and 
u is shorthand for cosO. dl/dz is the cosmological line 
element, Equation [S] 

The proper photon number density of a passively 
evolving component increases as (2 -1-1)'', and the gamma- 
ray attenuation per unit proper distance at a given en- 
ergy therefore increases by the same ratio. We can place 
the strongest constraints on reionization-era fields using 
the highest redshift gamma-ray sources available, pro- 
vided that these gamma rays satisfy the minimal energy 
condition discussed above. 



2.3 Spectral model 

Several emission processes can contribute to the cos- 
mological emissivity at the time of reionization. The 
fact that we are interested only in the broadband fea- 
tures of the spectrum, with integration over a range 
in redshift, will allow us to make some simplifying as- 
sumptions. A simple description of the stellar spec- 
trum at non-ionizing wavelengths as a thermal black- 
body is sufficient for massive pop-Ill stars that emit 
a large portion of their energy at ioniz ing wavelengths 
jSantos. Bromm fc Kamionkowskill2002l ). Ionizing radia- 
tion will undergo processing by neutral hydrogen within 
the host halo and eventually emerge at longer wave- 
lengths, and the fraction of emission that leaves the host 



halo (denoted by the escape fraction /esc) will be simi- 
larly afi^ected by the neutral hydrogen residing the inter- 
galactic medium (IGM). 

In this work, we focus on the contribution of stars 
with masses from 5 to 500 Mq to the cosmological back- 
ground radiation. For this purpose, we create rest-frame 
SED templates that include the stellar contribution as 
a truncated blackbody, plus Lyman-a, free-bound, and 
two-photon emission from reprocessed ionizing radiation. 
The mean free path for ionizing radiation is assumed to 
be short at all wavelengths we consider, and our tem- 
plates have no emission above the Lyman limit. We ig- 
nore the contribution to the spectrum from free-free in- 
teractions. While important at longer wavelengths, these 
only account for a small fraction of the emission at opti- 
cal and UV wavelengths that are the main concern in this 
work. Likewise we ignore the contribution of other recom- 
bination lines, which have only a relatively small effect 
that occurs at longer wavelengths. While these templates 
are not intended to accurately represent fine spectral fea- 
tures, they do provide a reasonable estimate of the en- 
ergy released into the cosmological background photon 
fields by early stars at rest-frame wavelengths from the 
Lyman limit to the optical, and can be raised or lowered 
to describe the output from a given global star-formation 
history. 

The IMF of low- and zero- metallicity stars is highly 
uncertain. Therefore, we consider two different possible 
IMFs in this work. The first is a simple Salpeter-like 
power law function from 5 to 500 Mp), while the second 
is a much more top-heavy Larson l|Larso nl ll998l ) IMF, 
defined as 



dN 
dM 



Af"^(l-f Af/A/c 



(11) 



which has a Salpeter slope at high masses (> Mc) 
and a flat slope at lower mass es, reducing the numbe r 
of low-mass stars. Following [Fernandez et al.l l|2010l ). 
and to allow a comparison with them in the next 
section, we set Mc = 250 M©. Our investigation of two 
different possibilities is motivated by the suggestion 
from simulations that there is a critical metallicity 
below which the formation of low-mass stars is strongly 
inhibited, leading to an IMF that is conside rably more 
top-h e avy than later population II stars (|Tumlinsonl 
I2OO6I : iNormanI bOOSj ). Whether this leads exclusively 
to very high mass stars (e.g. > 100 M©), and/or to a 
significant departure from a Salpeter -like slope in the 
high mass IMF {dN/dM oc M-^^^^ ISalpeteJ 'l955l is 
dispu te d by various au t hors jBromm. Coppi fc Larsor] 



20021: iBromm fc Loebl 12004: iTan fc McKed 12004 



TumUnson. Venkatesan fc Shull|[200i 'l. 

For the stellar component we apply a black- 
body spectrum down to a wavelength 912 A. Descrip- 
tions of the spectral temperatures, lifetimes, and num- 
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ber of ionizi ng photo n s rele ased by pop-Ill stars are 



taken from Schaerei 



As noted here and in 
the collective emission 



[Fernandez fc Kom atsu' (2006 
properties of high-mass stars per stellar mass are highly 
insensitive to mass above a certain threshold, ^ 50 — 100 
Mq. Above this threshold, stars are well described by 
an effective temperature ~ 10^ K, and release ~ 10®^ 
hydrogen-ionizing photons per solar mass over a lifetime. 
In the regime between the 5 Mq cutoff and this mass 
scale a dependence on the IMF still exists. In the case of 
our Salpeter IMF, about 81 per cent of mass in a zero- 
age population exists in stars of less than 100 Mq, and 
about half in stars less than 20 Mq , necessitating consid- 
eration of this deviation from the high-mass properties. 
In the case of the Larson IMF, < 100 Mq stars contain 
only about one third of all stellar mass. Integrated prop- 
erties with this IMF are therefore closer to what would 
be expected of an extreme IMF, populated only by > 100 
Mq stars. 



The Lyman-a line is the most significant spectral 
feature expected to appear in the reionization-era 
background. While higher-order Lyman-series lines are 
suppressed by repeated scattering event in the IGM, 
the photons in the a line will eventually escape out 
of resonance due to r edshifting and propag ate freely 
through the universe (|Loeb fc Rvbickil Il999t ). This is 
true regardless of whether the conversion of ionizing 
radiation to Lyman-a occurs inside the host halo or 
outside in the IGM. The overall spectrum is there- 
fore not significantly impacted by the ionizing escape 
fraction, which determines this ratio. While varying 
the escape fraction can alter the shap e of the angular 
pertu rbations in the near-IR EBL l|Fernandez et al] 
|2010[ ). the broadband properties of the UV emis- 
sion are not affected by th e details of this process 
l|Fernandez fc Komatsul [20061 ) . A simple assumption to 
determine the line strength is that some fraction of 
the energy emitted at ionizing wavelengths (which is 
excluded from the stel lar spectrum) eventually emerges 
cis Ly ma n-g emission ([Santos. Bromm fc Kamionkowskil 
2003: [Fernandez fc Komatsul " [20061) . Following 
Santos. Bromm fc Kamionkowskil ([20021 ) ■ Eq. (27), 



we assume 0.75 Lyman-a photons for each H-ionizing 
photon, plus an additional 0.53 for each He-II ionizing 
photon. The latter occur at about half the rate of 
hydrogen ionizing photons, for a high-mass star. 



Free-bound and free-free continuum radiation ac- 
count for a sub-dominant component of the emission in 
the UV. We include a description of the free-bound con- 
tribution in our spectral mod el, as per the derivation in 
[Fernandez fc Komatsul l[2006l ). As in that work, we as- 



sume an average gas temperature of T = 2 x 10 K. 
..(^,.) = 1.05xlO--|g^^,(.)Ax (12) 



gfb{n) (t>2{T) e fcrerg s"^ Hz"^ Mpc" 



Here p*{z) is the star- formation rate density (SFRD) in 
units of Mq yr~^ Mpc~'^, and dQu/diMQ) is the average 
number of ionizing photons released by the creation of 1 
solar mass of stars for an assumed IMF. R — 13.6 eV, gfb 
is the free-bound Gaunt factor and 4'2(T) is the recom- 
bination coefficient to all states other than the ground 
state, as described in Spitzcr (197§), Table 5.2. Both of 
these are close to one for our purposes. The summation 
over energy levels is truncated at n=3 in our calculation, 
as higher levels mainly contribute to longer wavelengths. 

A certain number of recombination photons will de- 
scend to the ground state by tw o-photon emis sion rather 
than Lyman-a. As described in lSpitzerl (|l978l ). two pho- 
ton emission is required for de-excitation of an electron 
in the 2'^S state, where single photon emission is prohib- 
ited. The probability of this occurrence is a weak function 
of temperature, and we use the valu e of 1 — zlc = 0.36 
fisted in Table 9.1 of lSpitzeil ([l978[ ). The emitted spec- 
trum, in terms of photon number, necessarily forms a 
symmetrical peak around 1/2 the Lyman-a energy (2432 
A or V2-y = t'Lya/2); the majority of the two-photon fiux 
emerges at UV wavelengths. The fitting function for the 
spectrum of two-photon e mission used here is taken from 
[Brown fc Mathewj ([l97d ): 



£^(1^) oc ——[1.307 - 2.627 ' 



+ 2.563 



V — V2'^ 



2V2 



7 



51.' 



2l^27 
2U2j 



(13) 



The total emissivity, which is then the input e(!/, z) 
for Eq. [T] is the sum of the Lyman-a , stellar, free-bound 
scattering, and two-photon contributions. We assume an 
instantaneous approximation for the emission of radi- 
ation, i.e., all light from a star is emitted the instant 
the star is created. Since the lifetime of high mass stars 
is very short compared to cosmological timescales (6.2 
xlO'^ yr for a 5 Mq star, 1.04 xlO'^ yr for a 15 M© star; 
[SchaereJ200i ). this is reasonable approximation that will 
have little impact on our results. 



3 RESULTS 



In i]3.1l we begin by setting some general limits on the 
pop-Ill contribution to the present-day EBL. This will 
be done by considering the highest energy photons de- 
tected from two distant gamma ray sources, and making 
the simple requirement that a passively evolving photon 
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population not create an optical depth that is greater 
than 1. In i|3.2l we will address the how these limits on 
the EBL can be translated into bounds on pop-Ill star 
formation. 



3.1 Limits on the EBL flux from reionization 

We find that a simple gamma-ray optical thinness cri- 
terion can significantly limit the high-redshift contri- 
bution to the observed EBL. In Fig. [1] we show the 
limits placed on EBL contributions from high redshift, 
based on two different gamma-ray sources that have 
been rece ntly observed with Fermi LAT. The first, GRB 
080916C (|Abdo et all 120091 ') was a brig ht gamma-ray 
burst with a measured redshift of z=4.35 (|Greiner et al.l 
[2009), which was observed by Fermi LAT and GBM on 
September 16, 2008. The spectrum of the source at high 
energy was found to be continuous over nearly 6 orders 
of magnitude, and the highest energy event was a 13.6 
GeV photon, corresponding to a rest-frame energy of 73 
GeV. The other source is flat-spectrum radio quasar PKS 
1502-1-10 6, which was detec ted during a flaring event at 
z= 1.839 l|Abdo et al.l2010bl) . The spectrum of this source 
was found to be most satisfacto rily fit by a l o g-parab olic 
function peaking near 1 GeV. lAbdo et al.l l|2010al ) re- 
ported a highest energy photon of 48.9 GeV from this 
source, corresponding to a rest-frame energy of 139 GeV. 

The thick colored lines in the figure delineate the al- 
lowed present-day fiux as a function of wavelength that 
does not produce an optical depth greater than 1 for the 
highest energy photons from each of these sources. These 
contours were created by considering a large number of 
hypothetical r-EBL SEDs, and excluding those for which 
the optical depth exceeded unity. These results are sen- 
sitive to the nature of the SED types considered - for 
this exercise we have used blackbody spectra at a va- 
riety of temperatures and normalizations. Because the 
optical depth from pair-production interactions (Eq. llOp 
depends on the integrated spectrum falling within the al- 
lowed energy range, the actual bound is dependent on the 
assumed type of spectral feature; a more sharply-peaked 
SED would be less constraining. Our use of simple black- 
body spectra avoids the introduction of any parameters 
into the results at this point. A couple of sample black- 
body spectra are shown on this plot, each at a normal- 
ization that would produce optical depth unity for the 
highest energy photon detected from GRB 080916C. 

This plot also shows the r-EBL SED models de- 
scribed in H2.3l that we employ in the next section, assum- 
ing a star formation rate density (SFRD) of 0.2 M0 yr~^ 
Mpc~^ from z = 15 t o z = 6. We hav e also plotted the 
EBL prediction from iFernandez et al.l (|2010l ) , when us- 
ing a star-formation efficiency of /, = 0.01. Differences 
in the shape of the prominent peak created by Lyman- 
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— 1 1 — I — I — r- 



— 1 1 1 — I — I — i- 



PKS 1502+106 



GRB OBogiec 




10* 

A{Angstroms) 

Figure 1. Upper limits placed on the contribution to the 
r-EBL by the assumption of the optical thinness of the uni- 
verse to the highest energy photons seen by Fermi from 
GRB 080916C (red. lAbdo et alj I2OO9I ) and PKS 1502+106 
(blue, lAbdo et al.]|2010bl ). The solid and dotted green lines 
are the SEDs of the p-EBL produced by pop-Ill stars form- 
ing with a Larson and Salpeter IMF, respectively, and with 
SFRD of 0.2 Mq yr-l Mpc'^ terminating at Zr = 6. For 
comparison, the Larson sp ectrum with a s t ar-for mation ef- 
ficiency of /« = 0.01 from [Fernandez et al" I 1 I2OIOI ) is shown 
as well (dashed black line). The two thin grey curves are 
examples of the thermal spectra used to create the opti- 
cal thinness r-EBL bound; here they are shown tangent to 
the limit for GRB 080916C. The thick blac k line is the pre- 
dicted p-EBL from the fiducial model of ISomerville et all 
l|201lh . Other grey points are total EBL measurements, 
including the lower l imits from number co unts from 
iMadau fc Pozzettil 1L2ODO") and lKeenan et alJ l|2010l ) (solid and 
open upward-pointing arrows), and sky photome t ry me asure- 
ments of DIRBE with sky sub traction: IWrightl 1I2OOII. solid 
square s) , ICambresv et all 1I2OOII. crosses), Levenson Sz Wrighli 
||2008| . solid di amond), iGoriian Wright fc Charvl (120001 
open hexagons ). IWright fc Reesel 1120001 open squares), and 



iLevenson. Wright fc Johnsonl | |2007| ' asterisks) . The small 
stars are direct m easurements from the IRTS experiment 
llMatsumoto et all [20051) . Some points are slightly offset to 
improve readability. 



a emission are due to the fact that star-formation is not 
constant in their model, and decreases with redshift. 

Our results here demonstrate the power of high red- 
shift gamma-ray source observations in limiting fiux con- 
tributions from high redshift. Both results constrain the 
flux at A ^ 1.4 fim from reionization to be a fraction of 
that produced below the redshifts of the two sources. 
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10 100 



E,(GeV) 

Figure 2. Attenuation for gamma-rays from tiigli redsliift 
sources as a function of observed energy. Solid black lines show 
the attenuation (e~'^) to gamma rays from sources at 2 = 2, 
4, and 6 (upper-right to lower left) predicted by the fiducial 
model of GSPDll (p-EBL; solid black line in Fig.[lJ. Dashed 
green curves are the attenuation for sources at these same red- 
shifts including the r-EBL component from Fig.[T] originating 
from high redshift pop-Ill stars with a Larson IMF and SFRD 
0.2 M0 yr~^ Mpc~^, ending at Zr = 6. 

Fig. [2] demonstrates the large increase in optical depth 
that the r-EBL components can create at high redshift 
compared to the p-EBL predictions in a typical model. 
By redshift 4, the Larson pop-Ill contribution shown in 
Fig. [T] produces a spectral cutoff at substantially lower 
energy than the GSPDll p-EBL alone, despite the fact 
that the latter dominates the former at all wavelengths 
in the z = EBL. 

3.2 High-redshift star formation 

By linking gamma-ray opacity to a particular high- 
redshift emission scenario, as described in ^2.31 we can 
constrain the amount of star formation consistent with a 
given set of gamma-ray observations. The limits derived 
will necessarily be dependent on the IMF applied, and 
the redshift at which zero-metallicity star-formation is 
assumed to terminate. 

3.2.1 The highest- energy photon method 

In this section, we describe how spectra from high- 
redshift gamma-ray sources observed by the Fermi LAT 



can be used to set limits on the background placed by 
pop-Ill stars. Because the most constraining sources will 
be those at the highest redshifts, we are interested in the 
small number of AGN and GRBs that have been detected 
at z 2 by Fer mi. A list of ' most c onstraining' sources 
was provided in lAbdo et all (l2010al ). To determine if a 
particular emission scenario is allowed or prohibited, we 
apply the highest energy photon method described in 
that paper. In this method, a power-law dN/dE = A 
is fit to the lower-energy signal from the source, which is 
assumed to be unaffected by photon-photon interactions 
(for PKS 1502-1-106, a log-normal parametrization was 
used). Once the parameters A and F are determined, the 
opacity induced by the background photon field, together 
with instrument parameters, are used to determine the 
probability distribution for the highest energy photon de- 
tected. If the likelihood of detecting this photon at or 
above the energy of the observed highest energy photon 
(Ehigh) is sufficiently small, then the opacity model can 
be excluded with a defined confidence level. This method 
provides a more robust exclusion of background models 
than the alternative of simply using a predefined value in 
gamma-ray opacity r at Ehigh, as we have done in Fig.[l] 
because it takes into consideration the normalization of 
the spectrum and the predicted number of observable 
gamma rays at and above Ehigh in the absence of EBL 
opacity, which may be much smaller or larger than unity. 

The expected number of photon counts above energy 
Ehigh is 

C{E)E-^e-^^''^dE (14) 

-^high 

Here A and F are the fit constants as previously de- 
scribed, while t(E) is the gamma-ray optical depth for 
the source at observed energy E, and C{E) is an energy- 
dependent factor describing the exposure of the source 
with Fermi LAT, i.e., instrument effective area integrated 
over time. The probability of observing zero counts is 
then described from Poisson statistics as 

P(0) = e-'^-'-^'-'f'-). (15) 

We have analyzed 5 high-redshift Fermi LAT sources 
selected from the 'most constraining' list of lAbd oetaLl 
iioioi), using publicly available data and software tools 
from the Fermi Science Support CenterQ The properties 
of the sources are summarized in Table [T] Our analy- 
sis uses the 'gtltcube' routine to determine source expo- 
sure, and the 'gtpsf function to calculate the energy- 
dependent point-spread function (PSF). We use the 
P6_V3 instrument response function, which is the default 
at the time of writing. Integrated data from the full in- 
strument lifetime is used for the AGN sources, and from 

^ http:/ /fermi. gsfc.nasa.gov/ssc/ 
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To - 5 to To + 135 s for GRB 090816C; To being the trig- 
ger time of the burst in the GBM instrument. Note that 
135 s is nearly 3 times the interval of GeV emission from 
the GRB presented in lAbdo et al] |2TOi). 

The photons within the 95 per cent containment re- 
gion of the PSF around the AGN source in each energy 
bin are then assumed to be associated with the source. 
For GRB 080916C, the 99 per cent containment region 
is used. To determine the normalization parameter {A) 
in Equation 1141 above, we use counts between energies of 
500 MeV and 5 GeV for GRB 080916C, and 1 and 10 
GeV for the four AGN. 



3.2.2 Combined probabilities method 

An alternative method is to consider the combined prob- 
abilities of having > 1 photons at or above the highest 
observed photon energy for each source, to put a stronger 
constraint on background fields. That is, we consider the 
combined probability 



PtotO 1)^Y[P^{> 1) 



(16) 



of detecting photons from all sources considered, and ex- 
clude background scenarios that lead to a Ptot less than 
a given value. However, a potential pitfall of this method 
is that Ptot can be considerably less than 1 even in the 
absence of a pop-Ill contribution; due to the impact of 
the p-EBL, which is quite uncertain at high redshifts. 
As discussed in G09 and GSPDll, this uncertainty in 
the UV background is a factor of several at z ^ 2. For 
high determinations of the UV background, like the fidu- 
cial model in GSPDll, the optical depth for the highest 
energy photons can be as high as 1, as shown in Table 
[T] The p-EBL effect is therefore not something that can 
be ignored here. Result derived using Eq. [TB] above will 
be necessarily be a strong function of the assumed star- 
formation model. In Table [1] we also show the expected 
number of photons expected above E^igh for each source, 
and the associated probability of detecting at least one 
photon at or above the highest photon energy. 

As we are interested in isolating the impact of the 
r-EBL, we first compensate for the impact of the p-EBL 
on our results by renormalizing the combined probability 
Ptot to 1, after convolving the observed spectrum with a 
given p-EBL opacity. Thus we define 

n,p(> 1) 



Ptot(> 1) = 



(17) 



where Pi^pEBL^^ 1) is the probability of 1 or more pho- 
tons from a source, after considering only an assumed 
p-EBL model, while Pi(^ 1) considers both background 
components. This change effectively isolates the impact 
of the r-EBL on photon detection probability. The result- 
ing limits on the reionization-era background are then 



only a weak function of the p-EBL; we find a variance of 
only about 10 per cent in our results for the pop-Ill star 
formation rate limits if the p-EBL flux is increased or 
decreased by a factor of 2. The following results assume 
the p-EBL predicted by the fiducial model of GSPDll. 

In Figs. [3] and |4l we show the overall limits on 
SFRD obtained using both the combined model, and 
a singular analysis using only GRB 080916C. Singular 
analyses with the other sources always lead to compara- 
tively weak upper bounds, 1 yr~^ Mpc~^ at the 
2(T level. A collection of observational data, as well as 
the model results for SGPDll are shown as well; how- 
ever most of these measurements reflect the luminosity 
of the brightest sources seen at these epochs, and are 
generally assumed to be due to population I/II stars 
in environments too metal-rich to support pop-Ill star 
formation. The measurement based on observed GRB 
rates are free from this bias, howev er they suffer from 
their own systematic uncerta inties (|Kistler et al.l l2009l : 
iBeckman fc Giammancoll2010l 'l. The pop-Ill SFRD in our 
model is assumed to be constant in time from redshift 15 
to the indicated cutoff redshift Zr. The choice of a func- 
tional form for the SFRD history is necessary and ar- 
bitrary here, and we have chosen a constant history for 
simplicity. In practice, the strongest effect on gamma-ray 
opacity will come from photons emitted near Zr, as evi- 
denced by the rapid increase in SFRD upper limits with 
increasing Zr in these figures. Results in these figures can 
therefore be considered as approximate SFRD limits at 
Zr. These limits are necessarily conservative, due to the 
renormalization that we have made in Eq. [T7|that limits 
the influence of the poorly-constrained p-EBL on our re- 
sults. If the p-EBL contribution to the UV is high, then 
the actual limits on pop-Ill stars after considering the 
total background photon population would be stronger 
than what is depicted in the figures. 

We have also shown in Fig. [3] the limits on pop-Ill 
growth that exist due to constraints on global metallic- 
ity and baryons available in collapsed structures. Pair- 
instability supernovae (SNpi) are responsible for a large 
release of metals into the IGM and ISM, which can 
potentially lead to confiict with the measured enrich- 
ment in high-redsh ift Lyman-a-forest observations (e.g. 
ISchave et al]|2003D . As a very simple enrichment model, 
we assume that all massive stars with 140 < M < 260 
Mq undergo SNpi upon deat h, and release half their mass 
into the IGM as metals (|Portinari. Chiosi fc BressanI 
Il998h . Enrichment from lower mass stars is ignored. For 
our Larson IMF, about 9 per cent of all stellar mass 
falls in the SNpi range, leading to a conversion of 4.5 
per cent of stellar mass into IGM metals with a sin- 
gle generation of stars. In Fig. [3] we show the level of 
global star-formation that between z = 15 and z = Zr 
would enrich the IGM to one-tenth solar metallicity. 
This is already considerably higher than measurements 
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Source ID 


Rcdshift 


Spectral Index 


^ norm 


Eh,ah (GeV) 


Tfid 




PpEBLO 1) 


GRB 080916C 


4.35 


-2.15 


36 


13.2 


0.34 


0.35 


0.295 


PKS0227-369 


2.11 


-2.56 


146 


31.9 


0.66 


0.20 


0.18 


PKS1502+106 


1.84 


-2.36 


2223 


49.2 


1.07 


1.06 


0.65 


PKS0805-07 


1.84 


-2.09 


973 


46.7 


0.96 


2.4 


0.91 


J1016+0513 


1.71 


-2.27 


322 


46.7 


0.83 


0.52 


0.41 



Table 1. Statistics for the five gamma-ray sources we consider in constraining pop-Ill star formation. Columns show the source 
name, redshift, and spectral index. Nnorm refers to the number of photons falling between 1 and 10 GeV (0.5 and 5 for GRB 
0809156C) that were used to determine the normalization of the source spectrum. -Ehigh is the highest energy photon received 
from the source, and Tfid is the optical depth of that photon produced by the p-EBL predicted in GSPDll. Nxjid{> E^igh) is 
the expected number of photon counts at and above E/^ig/^ after applying the optical depths from the GSPDll post-rcionization 
EBL, but before considering any contribution from pop-Ill stars, and PpEBhil? 1) is the probability of detecting at least one 
photon above Ef^igf^. 



such as those of Aguir re. Scha ve fc Theunj (|2002l ') and 
iBouche et aLl(|2007l ). which suggest a median enrichment 
of < 10"^ Zq at redshift 2, so our bound is quite con- 
servative. There are several ways in which this metal 
production could be suppressed. If the IMF were mod- 
ified from the Larson function assumed here, so as to 
reduce the number of stars forming in the SNpi mass 
range, metal production would be greatly diminished. 
Stars at masses higher than ^ 260Mq are generally be- 
lieved to follow a different evolutionary path than lighter 
stars and end with most of their metals inside remnant 
black holes (^Hc gcr fc Woos ley 2002 ). Shifting high- mass 
star production to either above or immediately below the 
SNpi mass range (i.e. 50 to 140 M©, or > 260 M©) would 
reduce IGM enrichment, while having little impact on the 
UV spectra produced per stellar mass. It is also possible 
that our assumption of all metals escaping to the IGM 
after an SNpi event is incorrect. With a Salpeter IMF, 
the mass fraction of baryons allotted to the SNpi mass 
range is much smaller, ~0.6 per cent, and the bound 
would thus be a factor of 15 higher. 



The number of baryons available at a given red- 
shift to form stars puts another constraint on any early 
star-formation model. The utilization of gas in collapsed 
structures can be quantified by the star-formation effi- 
ciency parameter /« . We take the total mass in all struc- 
tures above the molecular cooling mass scale as a function 
of redshift from fig. 1 of MS05, and use eq. (6) of that 
paper to compute the total mass conversion into pop-Ill 
stars for a given value /, ; 0.1 is used here. The correc- 
tion factor g in this formula is set to 0.65. While slightly 
higher masses of /, might be allowable, MS05 argue that 
/* ^ 0.3 is probably implausibly high. 



4 DISCUSSION AND CONCLUSIONS 

Our results indicate that the GeV sources seen by Fermi 
at 2 > 1.5 disfavor a scenario in which pop-Ill stars with 
a strongly top-heavy IMF are formed in copious num- 
bers (i.e. SFRD > 0.2 to 0.4 M© yr"^ Mpc"^) in the 
late stages of reionization, 6 < 2 < 8. At higher redshift, 
very high star-formation rate densities ( 5L 1) are disfa- 
vored by our result, and global limits on star- formation 
efficiency in proto-galaxies and limits on IGM metallic- 
ity impose a constraint at approximately the same level. 
Switching to a more moderate IMF, i.e. a Salpeter IMF 
with a cutoff at 5 M© , raises our limits by about a factor 
of 1.7. 

The results in Figs. [3] and |4] can be compared to 
those calculated in RKM09 for zero-metallicity stars 
that are based on low-redshift blazar limits on the lo- 
cal background (their fig. 9; the model with a = 10 and 
P = Q most closely resembles our simple step function ap- 
plied for our pop-Ill SFRD function). We place slightly 
stronger constraints on the SFRD at redshifts 7 and 8 
than this work. However, the significance of our strongest 
claims are only marginal, due to the very limited number 
of high-redshift photons that are relevant to our calcu- 
lation. The RKM09 result is based on an upper limit to 
the local background light of 5 nW m~^ sr~^ at 2 /im. 
As shown in Fig. [Jl our results most strongly constrain 
the EBL below ~ 1 /xm, and do not bound the EBL 
at co-moving wavelengths longer than 2 /xm. For simi- 
lar reasons, high redshift gamma-ray observations have 
little hope of providing meaningful constraints on a hy- 
pothetical contributio n to the IR background from dark - 
matter burning stars (|Freese. Spolvar fc Aguirr3l2008l l. 
Due to the relatively cool temperature of dark stars and 
the high redshifts at which they are theorized to exist, 
even an enorm ous contribution such as that in the max- 
imal model of iMaurer et all (|2010l ) would be at wave- 
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Figure 3. Upper limits placed on pop-III star formation rate 
density by our analysis of the spectra of high-redshift Fermi 
LAT detections. These results are for spectral output from 
stars assuming a Larson IMF. Solid red, blue, and green lines 
(bottom to top on right-hand side of plot) are limits set on 
SFRD for significance values of 2a, 3cr, and 5cr for the com- 
bined analysis (Eq. [T7J of all 5 sources listed in Table[T] These 
correspond to probabilities of 0.046, 0.0027, and 6.0 X 10~^ of 
finding the observed pattern of highest energy photons across 
all sources. The long-dashed orange line shows the 2a bound 
for the analysis of GRB 080916C, when considered alone. 
The dotted grey line at z ^ 6 shows the star formation rate 
level that would enrich the IGM to 0.1 Zq by Zr, while the 
dashed grey line corresponds to a star formation efficiency 
of /* =0.1 in collapsed structures; see text for more details 
on these criteria. We have also shown SFRD measurements 
at a variety of redshift. Because these measurements presum- 
ably include cooler pop-II stars at masses less than our cutoff 
of 5 Mq, they are not directly subject to our proposed lim- 
its on pop-III stars, but are included here for comparison. 
The solid black line is the s tar-formation his t ory pr edicted in 
the semi-analytic model of ISomerville et alj l|201 l|l: th is has 
been converted from a Chabrier IMF l lChabriedl20M) to a 
Salpeter with a multiplicative factor of 1.6, which is typical 
of the conversions described in lWilkins. Trentham fc Hopkina 
1I2OO8I '). Grw error bars without sym bols are from the com- 
pilation of iHopkins fc BeacomI |200^. Others measurements 
based on redshifted op tical/UV ligh t inclu de iBouwens et al.l 
1I2OO7I') (ope n circl es), IVerma et alj ||2007| ) (solid squares), 
Ouchi et al.l ll2p04|) (o pen squares), the upper limit from 



^^^^^c^_et_alj~ 2007 1 (downward-pointing triangle), and 



Thompson et all 1 200C ) (crosses) . The latter of these has been 



converted from a Chabrier IMF to Salpeter in the same man- 
ner as the semi-analytic model. Star-lik e symbols are results 
from GRB rate an alyses: IWang fc Pail f2009) (4-stars) and 
lYiiksel et al] |(2003) (6-stars). Some points have been shifted 
slightly for readability. 



Figure 4. As in the previous plot, but showing results for a 
Salpeter IMF (dN/dM oc M~^-'^^). Limits on star formation 
from IGM enrichment due to pair-instability supernovae are 
more than an order of magnitude higher in this case, and are 
not shown here. 



lengths too long for constraints to be derived, though 
this scenario would likely be difficult to reconcile with 
TeV-scale observations of lower-redshift blazars. 

We also strongly di sfavor a scenario in which a near- 
IR excess proposed by iMatsumoto et al.l (|2005l ) and a 
high EBL i nterpretation of the D IRBE 1.25 fim measure- 
ment (e.g. ICambresy et al.ll200ll ) are produced by early 
stars. Our results in Fig.[T]show that even a subdominant 
contribution to the EBL from high redshift at observed- 
frame wavelengths ^ 1 /im can be in disagreement with 
high energy data. We find good agreement with the result 
of ThoniDson ^et all (j2007l ). which disputes the detection 
of a high near-IR background, and attributes a possible 
contribution from high-redshift objects of 1-2 nW m~^ 
sr~^ at 1.4 to 1.8 /im using fluctuation analysis methods. 

In an analy sis of TeV blazars, 

lOrr. Krennrich fc Dwekl (|201ll ) recently found that 
the ratio of near- to mid-IR EBL flux needed to be 
significantly larger than prop osed in most recent mo dels, 
including those describ ed in iGilmore et al.l (|201ll ) and 
iDommguez et ahl (|201ll ). This work proposed a flux at 
1.6 nm of 17 ± 3 nW m~^ sr~^, which is above the 



level p roduced by resolved galaxies is lMadau fc Pozzett: 
200c ) and at the upper la bound of iKeenan et al 



20101 ). Our results do not strongly limit the possibility 



that an excess of a few nW m sr from high redshift 
could exist at or above this wavelength, but do disfavor 
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Figure 5. Limits on high-redshift SFRD that could be possi- 
ble with future gamma-ray observations. Contour lines on the 
plot show the redshift and highest photon energy of a gamma- 
ray source that would exclude a SFRD above the level indi- 
cated at redshift 6 with 2(t confidence. We make the assump- 
tion here that the gamma-ray source would have a spectral 
index F = —2.25 with no intrinsic curvature, and would have 
a normalization at lower energies such that one photon is ex- 
pected at or above the highest observed energy in the absence 
of background fields. Solid lines are results on SFRD for a 
Salpeter IMF, dotted lines are for a Larson IMF, and lines of 
each type are for limits of 1.0, 0.3, and 0.1 Mq yr~^ Mpc~^, 
from bottom left to upper right. Symbols on the plot indicate 
several of the most constraining sources observed thus far with 
the Fermi satellite. In addition to the 5 sources in Table [T] 
PKS 1144-379 and GRB 090902B are also shown. 



a scenario in which such an excess continues below a 
wavelength of ~ 1.2 /im. 

Our intent with this work is not just to compute the 
current limits on high-redshift star-formation available 
from gamma-ray observations, but also to demonstrate 
the usefulness of this method for future observations. 
Figs. [5] and [6] illustrate the SFRD limits that could be 
derived from future detections of high redshift sources 
with Fermi LAT or future telescopes. In these plots, the 
axes refer to the redshift and highest observed photon en- 
ergy of a hypothetical gamma-ray source. The source 
is then assumed to have a normalization at lower energy 
such that the expected number of photon counts at and 
above E^i is 1 [Nj:{> Ehigh) = 1] in the absence of any 
background field. The spectrum of the source is set here 
to -2.25, near the mean of the sources in Table[T] and the 
p-EBL is ignored. Given these parameters, the contours 



Figure 6. As in the previous figure, but for a cutoff redshift 

Zr = 9. 




Figure 7. Plot of the upper bounds on SFRD in two possible 
scenarios with future Fermi gamma-ray bursts, in the Larson 
IMF case. The solid lines show the limits from a GRB with the 
same redshift and spectral characteristics of GRB 080916C, 
but with a highest energy observed photon of 30 GeV (160 
GeV as emitted) instead of 13.2 GeV, in combination with 
the sources previously discussed. The dotted lines show a case 
with a GRB at 2 = 7, and a highest energy observed photon 
at 15 GeV (120 GeV emitted). Line colors arc as in Fig. [S] 
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on the plots show the source redshift and E-^ that would 
be required to place a given SFRD limit on pop-Ill star- 
formation at redshifts Zr = 6 and 9, with 2a significance. 
These contours are for limits derived based on a single 
source; combined limits for multiple sources like those in 
section 13.2.21 if available, would be somewhat stronger. 
In Fig. [T] limits based on two hypothetical high-redshift 
gamma-ray bursts are combined with the other sources of 
Table [T] This plot shows that new GeV sources, either at 
higher redshift than GRB 0809 16C, or at a similar red- 
shift with higher energy emission, could strongly limit a 
pop-Ill contribution to star-formation in the late reion- 
ization period. 

If the Fermi satellite remains in operation for its 
stated lifetime goal of ten years from its launch date, 
then its mission is currently less than one-third com- 
plete, and we can reasonably hope to see new GRB 
events or high-energy AGN photons that will strengthen 
our results. The Cherenko v Telescope Array (CTA; 
iThe CTA Consortium! lioiol ) is another possible source 
of constraining events. CTA will have a lower thresh- 
old energy than current-generation ground-based instru- 
ments, and may be able to detect sources at much higher 
redshift than currently achieved from the ground. Detec- 
tions with either of these instruments could potentially 
shed new light on star-formation in the reionization era. 
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